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Abstract 

The Douglas-Kroll transformed Dirac-Coulomb Hamiltonian is used to de- 
scribe scalar-relativistic effects in solids. A Hartree~Fock approximation with 
periodic boundary conditions makes it feasible to use methods originally de- 
veloped for atoms and molecules to solve the corresponding equations for 
crystalline systems. The implementation is realized within the CRYSTAL pro- 
gram. Scalar-relativistic effects in silver compounds of FCC structure such as 
AgF, AgCl, AgBr and Ag are investigated for the molar formation energy, the 
lattice parameter, the isothermal bulk modulus and the pressure derivative. 
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1. INTRODUCTION 



In the last decade, methods for treating relativistic effects in heavy atoms and molecules 
have been developed in order to extend the applicability of ab initio quantum chemical meth- 
ods to heavy atoms and molecules. Recently, different approaches and their realizations were 
compared as concerning the accuracy and efficiency of the method [T|]. A detailed discussion 
of the accuracy of different relativistic approximations by the investigation of hydrogenic 
energy levels can be found in Ref. 0. The Douglas-Kroll (DK) transformed Dirac-Coulomb 
Hamiltonian |^,^ is considered to be an efficient and reliable approximation for calcula- 
tions on relativistic systems. The theoretical framework of this approximation is based on 
quantum electrodynamics, and as one main consequence this relativistic Hamiltonian is 
variationally stable and therefore has normalizable bound-state solutions. The application 
of external-field projection operators which project onto the positive-energy states of the 
Dirac single-particle Hamiltonian leads to a two-component formalism which is separable in 
a spin-free and a spin-dependent part. There has been an extensive discussion of relativistic 
effects in the framework of this approach for atoms and molecules such as silver and gold 
0J^, gold hydride and gold fluoride [§] using the Hartree-Fock (HF) approximation as 
well as various methods including electron correlation effects and spin-orbit effects |p. 

There is also a plenitude of methods for treating relativistic effects in solids. The main- 
stream approach in this field being density functional theory (DFT) in the local density 
approximation (LDA). All these methods are, in principle, based on a relativistic extension 
of the Hohenberg-Kohn theorem [|l^ . Most of these approaches rely on the four-component 
treatment of a spherical symmetric potential (see Refs. [11,12] and references therein), which 
can be relatively easily solved using appropriate methods from atomic theory, extended to ef- 



fective one-particle equations given by MacDonald and Vosko [|^. One possibility to elude 
the shape approximation for the crystal potential are relativistic effective-core potentials 
(ECP) which are adjusted to (quasi-) relativistic atomic calculations (see Ref. [14] and 
references therein) and have to be transferable to the corresponding crystalline situation 
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The extension of HF methods in sohd-state theory to the relativistic case |[T^ is less 
obvious, however in practice, the treatment of potentials of non-spherical symmetry (and 
the non-local exchange potential) requires an effort which is much larger than in the non- 
relativistic case. Consequently, it seems to be a reasonable idea to use the DK transformation 
for the one-electron operators also in HF calculations with periodic boundary conditions. In 
this Letter we describe a corresponding implementation of the method in the CRYSTAL88 



program 19,20 and report first applications of the code. 



2. METHOD 

The spin-free part of the DK Hamiltonian is used to treat scalar-relativistic effects in 
solids. Using atomic units (with the electron mass m = 1), the relativistic operator has the 
following form: 

V:^{i) = - A [ Kxt(^) + ^ext(^) Ri]A-l {{E,, IVf (z)}, Wf{^}, (1) 



with A = R. 



2Ei ' mc'^ + Ei' 

The brackets {,} denote the anticommutator {a,b} = ab + ha. The operator Wf is an 
integral operator with kernel 

<(p.,p:)=a(r.-r:)a'^-^(p^'P^) 



Within this approach the use of the non-relativistic Coulomb repulsion instead of the trans- 
formed two-electron term is a good approximation if scalar-relativistic effects are considered 



An all-electron ab initio HF approximation for crystalline systems has been realized in the 
CRYSTAL88 program [^,^]. To construct a Fock operator in an atomic orbital (AO) basis 
(direct space representation) within this approach including the above relativistic operators, 
three main steps are necessary: 

1. In order to represent the operator RjV^xt(^)Ri in Eq- (0) in an AO basis, we have to 
calculate matrix elements of the operator pj Vext{i) Pi- 

2. The Bloch function (BF) representation of the kinetic energy, the external potential 
and the additional matrix elements mentioned above is used to get a representation of the 
relativistic operators in the BF basis. 

3. To get a representation of the relativistic operators in direct space, those quantities are 
Fourier transformed back to the AO basis. 



We make use of the fact that scalar-relativistic effects are short-range |]T6[. Therefore, the 
long-range tail of the electrostatic potential calculated in the CRYSTAL88 program by means 
of Ewald sums ||21]] is not relativistically corrected. Rather, the relativistic contribution of 
the external potential described by Eq. (0) to the Fock operator is restricted to a finite 
number of direct lattice vectors h within the spatial region that is treated by the direct 



lattice series of the Ewald potential function |^ - hereafter, dubbed "direct space potential 



zone" . 



Each AO is expanded in real spherical Gaussian-type functions (SGTF) p0| , |22| . We 
denote by [uj g] = Xn(r — g — s^) the uth AO of the primitive unit cell identified by the vector 
g located on an atom whose image in the primitive unit cell is assigned by the fractional 
vector Si^; g, h, ... label direct lattice vectors and Q = (n, I, m) denotes the quantum numbers 
characterizing the AO. For simplicity we write 1,2 for ui, U2 in the following. To evaluate 
the integrals, we calculate additional matrix elements of V^xt('^) — |r — h — s^l^^ 

in a basis of SGTF of the form 

i?(12g) =<V10|K^,(r)|V2g>. 

These integrals result from the operator R, which is equal to the linear momentum p = — zV 



times a p^-dependent kinematical factor. In practice, we calculate these derivatives with 
Cartesian Gaussian-type functions (CGTF) and transform the corresponding integrals to 
SGTF afterwards ||22|. 

In order to adapt the basis set to the translational invariance of the crystal, the AOs are 
used to build BFs pU|,pB|. It is convenient to use this basis to get a matrix representation 



of non-local operators such as E, A and {E + mc^)^^ because by definition this is the 
basis in which the matrix of one-electron operators is symmetry blocked. When the BF 
representation of the kinetic energy operator T 

f(12k) = S{12ky^/^ T(12k) 5(12k)^^/2 
is expressed in its eigenbasis 

f (12k) = [/(12k) A(12k) U\l2k) 

— the diagonal matrix A(12k) contains the positive real eigenvalues of T(12k) and the uni- 
tary matrix [/(12k) the corresponding eigenvectors — the BF representation of an analytic 
function /(T(12k)) corresponding to operators like £'(12k), yl(12k) and {E + mc2)^^(12k) 
is given by 

/(f(12k)) = [/(12k) /(A(12k)) [/t(12k). 

Here we assume that this relationship holds approximately also in the finite matrix repre- 
sentation given by the AO basis set, from which the BFs are constructed. The contribution 
of those operators to the relativistic external potential V^g'^^(12k) given by Eq. (|l|) — re- 
stricted to the direct space potential zone — is then realized by simple matrix multiplication 
for every k point within the irreducible wedge of the first Brillouin zone (IBZ). 

The relativistic external potential V^g'^^(12k) transforms like a totally symmetric op- 
erator. To transform this operator back to the AO basis, the Fourier transformation of 

C sf 

T/fj' (12k) restricted to the IBZ (see Refs. [20,23] and references therein) generates a matrix 



nf/(12g) = / Kff'''(12k) exp(^ k ■ g) dk 

JiBZ 



from which the full potential matrix V^g''^^(12g) is generated by applying the appropriate 
rotation matrices pO| , ^ . The relativistic operator £'(12k) is treated likewise. 

The Fock operator expressed in an AO basis can formally be written as sum of a one- 
electron and a two-electron part 

F(12g) = J/(12g) + 5(12g). 

The contribution of Eq. (0) to the Fock operator is related to H{12g) and does not include 
any relativistic correction of the long-range behaviour of the electrostatic potential Z(12g) 
because it is restricted to the direct space potential zone. We thus include relativistic effects 
due to the external potential by adding the difference of the relativistic potential and the 
non-relativistic potential calculated in this zone to the total non-relativistic potential as 
calculated in the CRYSTAL88 program. 

HiUg) = EiUg) + K^'^'(12g) - Kxt(12g) + ^(12g). 



Thus, the delicate treatment of the Coulomb series ||2T[] remains the same. As usual the total 
HF energy per primitive unit cell is given by 

e=^E ^(12g) (2//(12g) + i3(12g)), 

^ l,2,g 



where P(12g) denotes the density matrix of a crystal P4| . 



3. TECHNICAL DETAILS 

Scalar-relativistic effects in several simple silver compounds of FCC structure such as 
AgF, AgCl, AgBr and Ag are investigated using this method with the reoptimized uncon- 
tracted basis sets for F, CI and Br of Dunning and for Ag of Martin In crystalline 
compounds low-exponent functions are likely to introduce linear dependencies. Therefore, 
diffuse exponents (< 0.04) of the basis set for Ag are discarded. Concerning the question 
whether HF solutions for conductors are appropriate in spite of the vanishing density of 
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states at the Fermi level and the related problem of finding appropriate HF solutions of the 
free electron gas, we refer to Ref. [24] and references therein. Basis sets and details of the 
atomic and the CRYSTAL88 HF calculations using H^^ in the relativistic case and resulting 
total HF energies are given in Tables | and 0. 

All calculations are realized close to program limitations and whereever possible param- 



eters for high-precision calculations ||20| are used. To realize an accurate representation of 
the relativistic operators, it is convenient to use relatively large uncontracted basis sets. In 
order to obtain a reliable transformation from indirect space (k space) to direct space, the 
number of k points used to realize the integration over the IBZ should also be large. Testing 
the method described above for T(12k) and V^^^(12k) by comparing its IBZ integration with 
the given integrals of T(12g) and V^^^(12g), it turned out that 29 k points of the Monkhorst 
net [0 within the IBZ are sufficient to obtain an absolute error for the integrals which 
is less than 10~^. Moreover, the change in the total HF energy per primitive unit cell for 
metallic silver is less than 0.9 x 10~^ au if we use 72 k points instead of 29 k points within 
the IBZ for the relativistic version of the program. 

The geometry optimization is realized with the restart option GEOM of the CRYSTAL88 
program Two different equations of state (EOS) deduced by Murnaghan and Birch, 
respectively, are used to fit the total energy as a function of the lattice parameter. 
Though nearly all the results are only weakly dependent on the applied EOS — except the 
pressure derivatives of Ag and AgBr — the accuracy of the fit and especially related results 
such as bulk moduli and pressure derivatives have to be investigated carefully (for details 
see Table m). 

A corresponding ECP calculation of AgCl is repeated (see Table ^) to ensure 
compatibility of all used parameters. Since small-core pseudopotentials for Ag are used, it 
is assumed that frozen-core effects are small and relativistic all-electron HF calculations are 
comparable with corresponding ECP calculations. 
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4. RESULTS 



Although relativistic effects are widely discussed in literature (see Ref. [2S] and refer- 
ences therein), ab initio calculations of silver halides with FCC structure including scalar- 



relativistic corrections are rare |30] and in most cases pseudopotential approaches are used 
to consider those effects (see Ref. and references therein). By way of contrast, due to the 
simple structure of diamagnetic metals relativistic calculations of metallic silver are numer- 
ous (see Refs. [16,17] and references therein). Nevertheless, scalar-relativistic effects without 
any substantial approximation of the crystal potential have not been investigated before. 

The results of the investigated FCC crystals are collected in Table |T|. Properties such 
as the lattice parameters and the related molar volumes show small effects — they decrease 
up to 2.0% respectively up to 5.8% (for Ag) as compared to the non-relativistic values. The 
corresponding decrease of the nearest-neighbor distances lies between 0.8 pm for AgF and 



6.4 pm for Ag — relativistic bond-length contractions of similar molecules are larger |28,31 



Looking at the molar HF binding energies, a remarkable bond destabilization of the halides 
is observed — the molar HF binding energies are higher up to 8.9% (for AgF) as compared 
to the non-relativistic values. This trend is in accord with the finding in molecular com- 
pounds that noble-metal halides with a halogen more electronegative than the metal itself 
experiences a relativistic destabilization of the bond, along with a bond length contraction 



3^ . For Ag we find a large bond stabilization of about 34.5%. Bulk moduli and pressure 



derivatives change up to 7.5% (for Ag) respectively up to 8.2% (for AgCl) when the EOS of 



Murnaghan p3 is used. Corresponding fits of total HF energies versus the lattice parameters 
of AgBr and Ag are shown in Fig. 1 and Fig. 2., respectively. 

Although scalar-relativistic HF calculations are not directly comparable with calculations 
using scalar-relativistic ionic pseudopotentials in the framework of the LDA [|1^ due to 
the missing treatment of correlation effects, our results suggest that relativistic effects in 
4d metals are overestimated when ionic pseudopotentials are used (for a comparison of 
absolute values see annotation in Table fT|). This could explain why non-relativistic values 
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of 4d metals are slightly closer to experiment than scalar-relativistic ones (see Ref. [^] and 
references therein). 

In Table ^ we compare crystal properties of AgCl at various levels of theory. All calcu- 



lations are founded upon a periodic HF approach realized in the CRYSTAL program |19,33 



Results for different a posteriori density functional corrections and experimental values are 
included . We find remarkable agreement between relativistic HF calculations using H'^^ 
and relativistic ECP calculations using pseudopotentials of Hay and Wadt (see Ref. [29]) 
concerning the bulk modulus and the lattice parameter. Slight differences are observed for 
the molar formation energy and the pressure derivative. Corresponding fits of total HF 
energies per primitive unit cell versus the lattice parameter of AgCl are shown in Fig. 3. 



5. CONCLUSIONS 

We have calculated scalar-relativistic effects of crystalline properties such as lattice pa- 
rameters, molar formation energies, bulk moduli and pressure derivatives of various silver 
compounds with FCC structure by means of an implementation based on the Douglas-KroU 
transformation within a periodic Hartree-Fock approach realized in the CRYSTAL program. 

Good agreement between a relativistic effective-core potential calculation and the one- 
component Douglas-Kroll approach is obtained for most of the properties investigated for 
AgCl. Though scalar-relativistic effects of silver compounds are expected to be small, re- 
markable changes of the Hartree-Fock binding energy and of the bulk modulus are observed. 
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TABLES 



TABLE L Basis sets used for relativistic*^ and non-relativistic atomic Hartree-Fock (HF) cal- 
culations for the ground state of Ag, F, CI and Br and their total HF energies 



Atom Outermost exponents'' 
s p d 


Method State E*"* 


Ag Original basis: (15sl0p8d) [26] — > Final basis: (13s9p8d) 
Exponents deleted: 0.04 (s) 0.015 (s) 0.025 (p) 


0.74 (0.75) 0.81 0.40 (0.39) 
0.06 0.29 0.11 (0.12) 


HF limit [35] 2s -5197.6985 
nonrel. HF -5197.0041 
rel. HF^ -5310.3549 


F Original basis: (9s5p) [25] 


1.21 0.94 
0.36 0.27 


HF limit [35] -99.4094 
nonrel. HF -99.3955 
rel. HF^ -99.4821 


CI Original basis: (lls7p) [25] 


0.53 0.64 
0.19 0.18 


HF limit [35] -459.4821 
nonrel. HF ^p -459.4382 
rel. HF^ ^p -460.8422 


Br Original basis: (14sllp5d) [25] 


0.43 0.47 3.851 
0.16 0.14 1.32 


HF limit [35] ^P -2572.4413 
nonrel. HF ^p -2572.3577 
rel. HF^ ^p -2604.1851 



* Calculations using . 

^ Only reoptimized exponents (au) are tabulated. Exponents which differ for the relativistic calcu- 
lation are given in parentheses. 



TABLE II. Basis sets used for relativistic^ and non-relativistic CRYSTAL88 Hartree-Fock (HF) 
calculations of the FCC crystals Ag, AgF, AgCl and AgBr and their total HF energies (£'*°*) with 
reference to the primitive unit cell 



XJdOlO Oct 


Outermost exponents'' 


FCn rrvsital MptVinrI 


A r**^!! T r*'\r 


J-J 


for 


s 


P 


d 




parameter'' 


1 TP \ 


A IT 

Ag 


Original basis: (15sl0p8d) 


> r mai uasis. ^^losypoQj [see 


TqKIo Tl 




Ag 


0.77 


0.805 


0.385 


Ag nonrel. HF 




CI 07 ni Q1 

— oiy / .uiyi 




0.10 


0.22 


0.13 


Ag rei. tir 


e Q 7 7 1 Q 
D o < < io 


cqi n Q7c;i 

— OoiU.o t Oi 


Ag 


Original basis: (ISslOpSd) 


— > l^mal basis: (IdsypodJ [see 


lable iJ 




J_ 


OrigiiiMl 


basis: (9sr)p' 


[25] 








Ag 


0.77 


0.78 


0.43 










0.08 


0.24 


0.14 


x\gr llolirGl. rlr 


P, Q P, P, ^'3 
O D D lO 


— OZyO. 00410 


v 

r 


1.25 


0.95 




i\gr rei. rlr 


R Q p. p. ^o 
D o D D io 


— 04uy.y <o4 




0.38 


0.28 










Ag 


Original basis: (15sl0p8d) 


— > Jrimal basis: (13s9p8dj [see 


lable IJ 




V. . J. 


OrigiiiM] 


liasis: (Us?] 


>) [25] 


i |-' ' 1 1 "1 ^1 r"^ ^l 1 ^ ' f 1 1 ^ / ! ■> 1 rl 1 

J_ J.i l.<r li 1 In J.r> . \ .1. _L n J | ^ .H L / 






A ^ 

Ag 


0.79 


0.78 


0.40 










0.11 


0.24 


0.12 


/vgL^i noiiroi. rlr 


D o D D io 


— OOOO.OUi41 


\ji 


0.54 


0.70 




i\g*^i rei. Hr 


P Q P P ^'3 
D o D D io 


— / t i.o44li 




0.20 


0.20 


0.14 








Ag 


Original basis: (15sl0p8d) 


— > Final basis: (13s9p8d) [see 


Table I] 




Br 


Original 


l)asis: (] Isl; 


_ 1)5(1) [25 


] 






Ag 


0.76 


0.76 


0.40 










0.12 


0.11 


0.12 


AgBr nonrel. HF 


6 8 7 7 13 


-7769.5003 


Br 


0.46 


0.57 


3.851 


AgBr rei. HF^ 


6 8 7 7 13 


-7914.6672 




0.16 


0.19 


1.30 









^ Calculations using H^^ . 

^ Only reoptimized exponents (au) are tabulated. 

The meaning of the accuracy parameters is described in Ref. [20] and are given in the order 
described in the user manual of Ref. [19]. All values are obtained using 29 k points for the IBZ 
integration. 



TABLE III. Scalar-relativistic effects in tfie FCC crystals AgF, AgCI, AgBr and Ag^ 

Property Method AgF AgCl AgBr Ag'' 

00'=''' (pm) rel. HF 517.5 593.7 620.7 445.6 

nonrel. HF 519.1 599.0 626.0 454.6 

Vb (cm^mol-i) rel. HF 20.8664 31.5077 36.0048 13.3214 

nonrel. HF 21.0605 32.3591 36.9350 14.1450 

Eq (kJmoI-i) rel. HF -371.2 -385.9 -333.4 -53.0 

nonrel. HF -407.5 -417.7 -362.8 -39.4 

Bq^ (kbar) rel. HF 459 (459) 306 (309) 273 (272) 298 (295) 

nonrel. HF 460 (461) 297 (299) 259 (259) 322 (323) 

^0^^ rel. HF 4.09 (4.12) 5.03 (5.06) 5.69 (5.45) 1.40 (1.06) 

nonrel. HF 4.26 (4.31) 4.65 (4.75) 5.63 (5.52) 1.40 (1.18) 



^ Lattice parameter ao, molar Volume Vq, molar Hartree-Fock (HF) binding energy Eq with ref- 
erence to the HF ground-state energies of the free atoms and to the primitive unit cell, bulk 
modulus Bq and pressure derivative B'q. The relativistic all-electron CRYSTAL88 HF calculation 
using H^^ (rel. HF) is compared with a corresponding non-relativistic all-electron CRYSTAL88 
HF calculation (nonrel. HF). 

^ Comparing non-relativistic and scalar-relativistic results of Ref. [17] scalar-relativistic effects lead 
to a bond length contraction (nearest-neighbor distance) of 5.7 pm (1.9%), to a bond stabilization 
of -31.8 kJmor^ (11.0%) and to an increase of the bulk modulus of 210 kbar (19.8%) for Ag. The 
molar formation energy in Ref. [17] is calculated with reference to the ground-state energy of the 
free silver atom and to the Wigner-Seitz cell volume which is identical to those of the primitive unit 
cell. The corresponding values for Ag in Table HI are: bond length contraction 6.4 pm (2.0%), bond 
stabilization —13.6 kJmol"^ (34.5%), decrease in the bulk modulus 24 kbar (7.5%) (percentage 
with reference to non-relativistic values). 

Parameters of the restart option GEOM of the CRYSTAL88 program are for AgF (5.00), AgCl 
(5.90), AgBr (6.20) and Ag (4.45). The accuracy parameters are the same as in Table H. 
^ Ranges of bond distances used for fitting the EOS: AgF [490 - 540], AgCl [560 - 630], AgBr [605 
- 670] and Ag [425 - 510] (all values in pm). The ranges have been chosen large enough to ensure 
numerical stability, and small enough to exclude inappropriately described regions of the potential 
curve at large bond distances, taking the inadequacy of the closed-shell HF approximation for near 
degenerate situations into account. Different EOSs [27] lead to same results for oq - otherwise 
values with an EOS of Birch [27] are given in parentheses. All mean square deviations of the fits 
are less than 0.6 x 10~^. 



TABLE IV. Crystal properties of FCC AgCl at various levels of theory^ 



Method E^°^ (a.u.) Eq (kJmori) oq (pm) Vq (cm^mol-i) Bq (kbar) Bq 



nonrel. HF^'^ 


-5656.6014 


-417.2 


599.0 


32.3591 


297 (299) 


4.65 (4.75) 


rel. HF'^'^ 


-5771.3441 


-385.7 


593.7 


31.5077 


306 (309) 


5.03 (5.06) 


rel. ECP^'^ 


-159.7086 


-376.5 


594.3 


31.6000 


307 (312) 


5.38 (5.44) 


rel. ECP^ 


-159.708 


-375.4 


594.0 


31.5532 


330 




DFTie 


-160.656 


-554.0 


570.0 


27.8810 


510 




DFT29 


-160.870 


-514.6 


561.0 


26.5811 


560 




exp.*^ 




-519.8 


551.0 


25.1848^ 


535 [417] 


[7.00] 



^ Total energy with reference to the primitive unit cell, molar formation energy Eq with 
reference to the ground-state energies of the free atoms and to the primitive unit cell, lattice 
parameter oq, molar volume Vq, isothermal bulk modulus Bq and pressure derivative B'q. 
^ Non-relativistic all-electron Hartree-Fock (HF) calculation using the CRYSTAL88 program. 

Relativistic all-electron Hartree-Fock (HF) calculation using the CRYSTAL88 program withiJ^^i. 

Relativistic effective-core potential (rel. ECP) calculation using the CRYSTAL92 program and 
ECPs for Ag and CI given in Ref. [33]. 

^ Relativistic effective-core potentials (rel. ECP) results calculated in Ref. [29] using nearly the 
same ECPs as in footnote d. (For Ag a (5s5p4d)/[3s3pld] GTF basis set is used instead of a 
(5s5p4d)/[3s3p2d] GTF basis set with ^ = 0.2108 for 4d instead of ^ = 0.2079 for 4d). 
^ See footnote d of Table HI. 

§ Correlation contribution calculated in Ref. [29] with a posteriori density functionals of Colle and 
Salvetti (DFTl) and Perdew (DFT2). 

^ Estimated experimental values (see Ref. [29] and references therein). The estimated values for 
£^0) and Bq are an extrapolation of experimental results to T=0 K including a zero-point energy 
correction. All values are zero-pressure values. A zero-pressure value of the pressure derivative Bq 
for the Bl phase of AgCl at T=0 K is not available. The estimated zero-pressure value Bq at room 
temperature and the corresponding pressure derivative Bq (see Ref. [34]) are given in parentheses. 
' Vq is deduced from oq. 



FIGURES 




FIG. 1. Total Hartree-Fock (HF) energy versus lattice parameter of FCC AgBr. The fits using 
an EOS given by Murnaglian [27] are shifted by 7769.0 au (nonrel.) and by 7914.165 au (rel.). Ab- 
breviations used: rel. (relativistic CRYSTAL88 HF calculation using nonrel. (non-relativistic 
CRYSTAL88 HF calculation). 




FIG. 2. Total Hartree-Fock (HF) energy versus lattice parameter of FCC Ag. The fits using an 
EOS given by Murnagtian [27] are shifted by 5197.0 au (nonrel.) and by 5310.35 au (rel.). Abbre- 
viations used: rel. (relativistic CRYSTAL88 HF calculation using H^^), nonrel. (non-relativistic 
CRYSTAL88 HF calculation). 
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FIG. 3. Total Hartree-Fock (HF) energy versus lattice parameter of FCC AgCl. The fits using 
an EOS given by Murnaghan [27] are shifted by 5656.0 au (nonrel.), by 5770.7 au (rel.) and by 
159.062 au (rel. ECP.). Abbreviations used: rel. (relativistic CRYSTAL88 HF calculation using 
H^^), nonrel. (non-relativistic CRYSTAL88 HF calculation), rel. ECP (relativistic effective-core 



potential CRYSTAL92 HF calculation [33]). 



